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A  general  purpose  dyadic  diffraction  coefficient  for  plane  and 
spherical  waves  obliquely  incident  on  a  perfectly-conducting  circular  cylinder 
is  derived.  This  diffraction  coefficient  is  then  employed  to  calculate 
the  radiation  pattern  of  a  LAMPS  antenna  with  a  cylindrical  mast 
positioned  in  front  of  it. 


DD  I  j^M,j  1473  edition  of  NOV  A  5  IS  OBSOLETE 


UNCLASSIFIED 
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I.  INTRODUCTION 

The  development  of  a  computer  code  which  can  be  used  to  calculate 
the  effects  of  a  circular  cylinder  on  the  radiation  pattern  of  an  aperture 
antenna  is  described  in  this  report.  The  motivation  for  this  research  is 
to  investigate  the  pattern  degradation  of  shipboard  antennas  caused  by 
nearby  masts. 

The  report  is  divided  into  two  parts.  The  first  part  is  concerned 
with  the  derivation  of  expressions  for  the  fields  diffracted  by  a  rela¬ 
tively  small  (in  terms  of  A)  circular  cylinder.  The  resulting  solutions 
for  this  diffracted  field  are  put  into  the  GTD  format  so  that  the  cylinder 
diffraction  coefficients  can  be  identified.  Thus,  the  results  can  be 
incorporated  into  a  general,  ray-optical  analysis  of  a  radiating  system. 

The  second  part  of  the  report  is  concerned  with  using  the  above  results 
to  predict  the  degradation  of  the  pattern  of  the  LAMPS  antenna  caused  by 
the  presence  of  a  nearby  cylinder.  The  LAMPS  antenna  is  a  34"  parabolic 
dish  with  a  nominal  freguency  of  4.6  GHz.  The  resulting  computer  code  is 
capable  of  calculating  the  radiation  pattern  of  the  LAMPS  antenna  both  with 
and  without  the  diffracting  cylinder  present,  so  that  the  effects  of  the 
cylinder  on  the  radiation  pattern  can  be  easily  seen.  While  the  present 
computer  code  is  limited  to  calculating  patterns  in  a  plane  perpendicular 
to  the  cylinder  axis,  with  the  cylinder  being  parallel  to  the  plane  of  the 
aperture,  the  theory  presented  herein  is  adeguate  to  provide  for  an  exten¬ 
sion  of  the  code  to  the  more  general  case. 

II.  CYLINDER  DIFFRACTION  COEFFICIENTS 

The  goal  of  this  derivation  is  the  development  of  diffraction  coeffi¬ 
cients  which  can  be  used  to  find  the  diffracted  fields  due  to  an  incident 
wave  interacting  with  a  perfectly-conducting  circular  cylinder.  Two  cases 
will  be  considered:  the  soft  boundary  case  (H  field  perpendicular  to  the 
cylinder  axis)  and  the  hard  boundary  case  (E  field  perpendicular  to  the 
cylinder  axis).  The  method  will  be  guite  straightforward;  a  canonical 
problem  will  be  solved  and  asymptotic  approximations  will  be  made  so  that 
the  resulting  solution  can  be  interpreted  ray-optically.  Then  this  solution 
will  be  compared  with  a  GTD  solution  to  obtain  the  unknown  diffraction 
coefficients . 

It  will  be  seen  that  the  resulting  diffraction  coefficients  are 
analogous  to  those  encountered  in  wedge  diffraction,  which  is  not  sur¬ 
prising  since  in  both  cases  one  is  dealinq  with  a  line  of  scattering  (the 
axis  of  the  wire  or  the  edge  of  the  wedge). 
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B.  Soft  Diffraction  Coefficient 


The  geometry  of  the  canonical  problem  we  wish  to  consider  is  shown 
in  Figure  1.  The  field  incident  on  a  perfectly-conducting  cylinder  of  radius 

z 


Figure  1.  z  directed  current  moment  (dipole  in  the  presence  of  a  circular 
cylinder  centered  on  the  z  axis. 


a  is  due  to  a  z-directed  electric  current  moment  pe=H  at  p=Pq,  t=0,  z=0.  TF 
current  I  extends  over  the  incremental  distance  £«A.  The  cylindrical  dif¬ 
fraction  point  Qc  is  located  on  the  axis  of  the  cylinder  at  the  point  where 
fc><H=6 i .  The  cylinder  (or  wire)  like  the  edge  of  a  wedge  is  a  line  of 


scattering,  so  the  angle  of  diffraction  b^  equals  the  angle  of  incidence  e  j 
The  field  point  is  located  using  the  ordinary  spherical  coordinates  (r,e,<j>) 

Considering  the  problem  ray-optical ly  [2],  we  can  express  the- 
component  of  the  diffracted  electric  field  as 


Ed(0  -  Ef(qc)0 sj 


S'  _ 


e 


-jks 


(1) 


in  which  E^Qc)  is  the  incident  electric  field  at  Qc,  Ds  is  the  soft 
diffraction  coefficient,  and  s'  is  the  caustic  distance,  see  Equation  6.42 
of  [2].  Note  that  E1  is  in  the  plane  of  incidence  formed  by  the  cylinder 
axis  and  the  ray  incident  from  the  current  moment  (or  dipole)  to  Qc.  The 
above  expression  is  analogous  to  that  used  to  calculate  the  fields  dif¬ 
fracted  by  a  wedge;  the  diffraction  coefficient  remains  to  be  found. 


When  Qc  is  in  the  far-zone  of  the  current  moment,  the  incident  field 


Ei(Qc)  -  -j 


kpen 

4tt 


sine^ 


(2) 


in  which  k  =  2tt/>  is  the  wave  number  and  n  is  the  characteristic  impedance 
of  free  space.  Combining  (1)  with  (2) 


Ea(s)  =  -j 


kPeh 


4tt 


-  sine^D 


•jk(s'+s) 


1  s  -Jss 1  ( s+s 1 ) 


(3) 


One  can  obtain  an  eigenfunction  solution  for  this  same  geometry. 
When  the  field  point  is  in  the  far  zone,  the  e  component  of  the  total 
electric  field  is  given  by 


kp  n  sine  -jkr  jkp  sinocost 
Eo  ■  J  lil  -  V 


where 


«  J  (ka  sin  o) 

SP  =  I  f-J  Jrr\ - Hy;(kP„sino)cos  m* , 


e 

m=u 


(ka  sine)  m 


the  phase  reference  is  the  origin,  and 

1 ,  m=0 

2,  m^O. 


m 


(4) 


(5) 


In  the  far-zone  note  that  e=e.=e,. 


The  result  given  in  (4)  and  (5),  which  is  derived  in  the  Appendix, 
has  been  obtained  previously  by  reciprocity  arguments  [3],  the  Hertz 
vector  [4]  and  the  Green's  function  method  [5]. 

In  (4)  it  is  apparent  that  the  first  term  is  the  incident  field  of 
the  dipole,  so  the  second  term  must  be  the  field  diffracted  from  the 
cylinder.  Referring  the  phase  of  the  second  term  to  Qc  with 

r  s  +  zQ  cosO  ,  (7) 


4tt 


si  no 


-jkz  cose  e-^3 
e  s 


(8) 


The  above  equation  may  be  compared  now  with  the  far-zone  form  of  (3)* 


Ed(s) 


J  4  IT 


sine  Ds 


(9) 


to  obtain 


-jkz  coso 
e  0 


(10) 


In  this  form  Ds  is  dependent  on  the  distance  s'  from  Qc  to  the  source, 
the  distance  z0  from  Qc  to  the  origin,  and  the  distance  P0  from  the 
origin  to  the  source;  where  P0  appears  in  the  expression  for  Se-  We 
will  next  simplify  the  above  equation  for  Ds  to  remove  the  dependence 
on  these  distances,  so  that  Ds  will  depend  only  on  the  cylinder  radius 
a,  a,  and  the  aspects  of  incidence  and  diffraction  as  viewed  from  Qc. 

A  case  commonly  encountered  in  practice  occurs  when  Pq  >■>  a,  i.e., 
the  dipole  is  many  radii  removed  from  the  cylinder.  For  this  situation 
we  can  introduce  the  large  argument  asymptotic  approximation  of  the 
Hankel  function, 


H^)(kposin0) 


1.  — 
4* 


Ji 


•m 


-jkp0sine 


kp0sine 


(11) 


into  the  equation  for  Se.  In  the  above  approximation  the  argument  must 
be  much  larger  than  the  order;  however  this  condition  is  met  for  the 
terms  in  the  series  which  contribute  significantly  to  Se-  Thus  Se 
reduces  to 


*t)ote  that  the  E-field  unit  vectors  defined  in  Figure  1  are  in  the 
-o  direction. 
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2j 


'f  Tkp  sino 


■ikp  sinb  “  ,  Jp(ka  sinfl) 

£  k  - 1  )  Tq  y  COS  1.1$ 

pfO  '  H'd'(ka  sin-0 


n 


U) 


rployirq  the  abov<:>  equation  in  (10)  and  utilizing  the  relationships 


s'  =  p  /sim  =  p  sine  +  z  cose, 
o  o  o 


(13) 


it  is  found  that 


i)  =  -  M  4 

s  V  tt  k 


71  k  sin9d  n=0 


“■  J  ( ka  sine  ) 

l  ‘m*-’) '  0(?ir~ —  cos 


Hv  ;(ka  sine  ,] 
m  d 


(14) 


when  pq  * •-  a.  This  is  the  desired  form  for  the  soft  diffraction  coeffi¬ 
cient,  since  it  depends  only  on  a/x  and  the  aspects  of  incidence  and 
diffraction.  It  may  be  used  in  (1)  even  though  the  field  point  is  not 
in  the  far  zone,  where  s  ■  s';  it  is  sufficient  that  Is  sinfe^  1. 

This  explains  why  has  been  used  in  place  of  fe  in  (14).  For  a  very 
thin  wire  where  ka  -  1,  the  formula  for  Ds  simplifies  further  to 

-J  i 

D  =  .-.e -  — - — - -  (15) 

s  \f 2t- k  sinQ^  in  2-in(v  ka  sined)-j  j 

in  which  v  =  1.781.  It  is  seen  that  Ds  is  independent  of  $  for  a  sufficiently 
thin  cy I i nder  (wi re) . 


It  is  evident  that  Ds  given  by  (14)  is  much  more  convenient  to 
use  than  the  form  given  earlier  in  (10).  However  in  cases  where  p0 
is  not  much  greater  than  a,  it  is  desirable  to  seek  a  correction  to 
(14)  which  introduces  PQ  in  a  relatively  simple  way. 


Let  us  consider  the  Debye  or  tannent  approximation  to  the  llankel 
f’/nction.  This  is  a  more  accurate  asymptotic  approximation  and  is  not 
limited  to  the  case  where  kPQsi ne-’-'n.  The  Debye  asymptotic  anproximati on 
is 


n  '  1 


-A. 


’Txsina 


e-jx(sinu--acosa) 


H6) 


where  x  cosa  =  m. 
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to  m) . 


The  following  approximations  are  valid  for  m/x  <  1  (but  x  not  close 


ttx  sina  v  TTX 


,  m  1  2 
’ + 1  <f  h 


x  sina  v-  x  -  ni 


(17a) 

(17b) 


x  a  cosa  %  m  -k  -  m  2 


(CL)  +  i  (CL) 

l2x'  3  2x 
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(17c) 


Substituting  the  above  into  (16)  and  noting  that 


jr.i  2 


.m 


=  J 


H*2)(x)  '  {%  J”  e-J* 


'  *\  <!> 


21 


♦  \  &> 


31 


(18) 


It  is  evident  that  the  term  in  the  bracket  and  the  following  exponential  term 
are  approximately  equal  to  unity  for  x->>m;  hence  the  term  preceding  the 
bracket  is  the  large  argument  approximation  (see  (11)  where  x=kP0sine).  Since 
the  series  representation  of  Se  given  in  (5)  converges  for  m%2kasin0,  it 
follows  from  the  condition  x>n  that  ^Q^a,  since  x=kP0sine  here.  The 
condition  P0>2a  is  usually  met  in  practice. 


The  relative  size  of  the  correction  terms  in  (18)  which  contain  the 
ratio  m/x  (with  m%2  ka  sine  and  x=kP0  sine)  indicates  whether  the  form  of 
Ds  with  the  large  argument  approximation  for  the  Hankel  function  is  valid. 
If  the  correction  terms  are  too  large  to  be  neglected  and  pn>2a»  the 
Debye  form  of  the  Hankel  function  in  (18)  can  be  used  in  (10)  with  the 
result  that 


rx-r  ,  J  (ka  sine  ,) 

n  I  2j  1  V  /  ,  m  _ 

s  =  *  Uk  sined  J0  em  "  ^(ka  sine 


cos  m<f>  x 


1  +1‘  m 


2  ^p0sined 


r 


exp<-jri 


m 


+  1 


-1 
1  i 


^2kpQsin0  J  3  \2kp„sin0 


osin9d  )_ 


(19) 
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B.  Hard  Diffraction  Coefficient 

One  can  proceed  in  a  similar  manner  to  find  the  hard  diffraction 
coefficient  Dh  by  replacing  the  z-directed,  electric  current  moment  by 
a  z-directed  magnetic  current  moment.  Note  that  E1  is  now  perpendicular 
to  the  aforementioned  plane  of  incidence.  For  this  case  the  .'-component, 
of  the  diffracted  field  is  given  by 


Ed(s)  =  E1(gc)Dh|il^f?TT  e 


in  which  E  (Qc)  is  the  (scalar)  incident  electric  field  at  the  point  o’ 
diffraction  Qc .  Following  the  same  steps  used  to  obtain  (10)  for  the 
soft  diffraction  coefficient,  we  find  that 


0 


h 


is~  e-jks' 


s 


m 


- jkzQcosO 
e 


(21 


where  the  expression  for  Srn  derived  in  the  Appendix  is 


m  J‘(ka  sin") 

j  jjy, - (kn  sino)cos  m<t 

H^>  (ka  sino)  m  0 


in  which  the  prime  denotes  differentiation  with  respect  to  the  argument. 

For  the  case  where  r0-’:’a,  the  large  argument  approximation  for 
Hr|)£-)(kn0sin'J)  may  be  used  as  before  to  obtain 


(U 

21  JL_  j 

irk  sine  L 


d  m=0 


J '  ( ka  si  no  ) 
,m  m  d 


Hu'  (ka  sino. 
m  d 


cos  mi 


which  compares  with  (14)  for  the  other  polarization  (soft  case).  For  a 
very  thin  wire  where  ka  ^  1,  the  above  expression  simplifies  further  to 


Comparing  the  above  equation  with  (15),  it  is  clear  that  iDh!""  Ds  1  ,  as 
one  would  expect  for  the  thin  wire  case.  Furthermore  one  notes  that 
there  is  null  in  the  diffraction  pattern  at  =  120’. 
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If  i'0  is  not  ""a  but  o02a,  we  may  employ  the  Debye  approximation 
of  (kP0sinO)  with  the  result  that 


d  =  -  /n  _ j ?  (., 

h  Jt Tk  sin1*,  Ln  1 


d  m=0 


J  (ka  sin*,) 
v  m  nr  d; 

)  -ryy f - COS  m  !  x 

Hm  (ka  sin  'd> 


which  compares  with  (19)  for  the  other  polarization. 

The  scalar  diffraction  coefficients  may  be  employed  in  a  dyadic 
diffraction  coefficient  of  the  same  form  as  that  derived  for  edge 
diffraction  [1].  As  mentioned  earlier,  the  plane  of  incidence  is  defined 
to  be  the  plane  containing  the  incident  ray  to  the  diffraction  point  Qc 
and  the  axis  of  the  cylinder,  and  in  addition  we  may  introduce  a  plane 
of  diffraction  defined  to  be  the  plane  containing  the  diffracted  ray 
and  the  axis  of  the  cylinder.  The  angle  between  these, two  planes  is  i>. 
Let  us  now  decompose  the  incident  electric  field  into  E 1  ( Qc )|(  and  E “*  ( Q c ) ,  , 
components  parallel  and  perpendicul  ar  to  the  plane  of  incidence;  simi¬ 
larly  let  us  separate  the  diffracted  electric  field  nto  Ed(s)„  and  Ed(s)x, 
components  parallel  and  perpendicular  to  the  plane  of  diffraction.  Then 
in  matrix  form 


where  the  directions  of  the  parallel  field  components  are  shown  in  Figure  1, 
and  the  direction  of  the  perpendicular  component  in  each  case  is  determined 
by  taking  the  vector  product  of  the  parallel  unit  vector  (see  Figure  1)  with 
a  unit  vector  in  the  direction  of  propagation.  We  conclude  this  section 
by  presenting  some  numerical  results  for  the  magnitude  of  Ds  and  Dh  as  a 
function  of  4>  for  various  cylinder  radii  from  ka  =  .25  (a  =  .04*)  to 
ka  =  4  (a  =  .64*);  0<j  =  90°  in  each  case.  Equations  (14)  and  (23)  were 
used  to  calculate  the  curves  shown  in  Figures  2  and  3.  The  number  of 
terms  included  in  the  series  varied  with  cylinder  size.  For  the  two 
smaller  cylinders  4  terms  were  used,  for  the  ka  =  1  cylinder  5  terms, 
and  for  the  ka  =  4  cylinder  8  terms.  In  line  with  the  earlier  discussion, 
larger  diameters  could  be  accommodated  at  the  cost  of  including  more  terms 
(approximately  2ka  sinOj)  in  the  series.  It  is  seen  that  the  curves  for 

ka  =  .25  approximate  the  behavior  predicted  by  (15)  and  (24)  which 
described  earlier. 
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III.  MODELING  AN  APERTURE  ANTENNA 


Let  u'j  now  proceed  to  the  computer  code  developed  to  calculate  tne 
cylindrical  aperture  blockage.  The  problem  which  can  be  handled  with  the 
present  program  is  shown  in  Figures  4  and  6.  The  equivalent  currents  in 
the  aperture  are  modeled  with  z-cnrected  electric  and  magnetic  dipoles 
situated  in  the  YZ-plane.  The  magnitudes  of  the  various  dipoles  can  be 
adjusted  so  as  to  model  tapered  or  other  non-uniform  aperture  distributions. 
The  cylindrical  scatterer  of  radius  a  is  located  at  p  -  po,  i>  -  io*  and  is 
parallel  to  the  z  axis.  The  program  is  capable  of  computing  the  far-field 
patterns  in  the  XY  (i.e.,  fl  =  90’)  plane,  both  with  and  without  the 
cylinder  present.  Thus  only  a  one-dimensional  model  of  the  aperture  is 
required. 


The  dipole  moments  Idl  and  Kdl  for  the  equivalent  electric  and  magnetic 
dipoles  are  obtained  from  the  electric  field  aperture  distribution  as  shown 
in  Figure  5.  The  aperture  is  divided  into  thin  vertical  strips  of  width  W, 
where  the  W  is  chosen  small  enough  so  that  the  relative  phase  of  the  aperture 
fields  are  essentially  constant  across  the  strip  for  all  pattern  aspects. 

From  the  uniqueness  theorems  we  know  that  only  the  y-  and  z-  (tangential) 
components  of  the  E  field  (or  H  field)  need  to  be  specified  in  the  aperture. 
For  a  planar  aperture,  assuming  that  the  fields  in  the  aperture  plane  out¬ 
side  the  aperture  region  are  negligibly  small,  we  have,  from  the  equiva¬ 
lence  principle, 


Jz(y,z)  =  -2  E2(y,z)/n  (27) 


Kz(y,z)  =  -2  Ey(y,z) 


(28) 


where  J  and  K  are  the  equivalent  electric  and  magnetic  surface  currents  in 
the  aperture,  the  subscripts  denote  the  y  and  z  components,  and  n  is  the 
free  space  characteristic  impedance. 


We  can  easily  obtain  the  electric  and  magnetic  dipole  moments  cor¬ 
responding  to  the  equivalent  currents  for  any  aperture  strip.  For  the 
dipoles  located  at  y=yn  we  have 


Idl 


f2b  ,VW/2 

JJy,z)ejkz  cos0  by  dz 

l-zb  Jyn-w/2 


(29) 


Kdl 


yn+w/2 

K  (y,z)e^kz  cos9  dy  dz 

yn-w/2 


(30) 
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Figure  4.  The  present  computer  program  is  capable  of  calculating 
the  far  field  pattern  as  a  function  of  $(0=90°)  for 
the  geometry  shown.  The  aperture  antenna  is  modeled 
by  the  z  directed  magnetic  and  electric  dipoles  located 
along  the  y  axis.  The  dipole  magnitudes  need  not  be 
identical,  so  that  tapered  aperture  distributions  may 
be  modeled. 


The  far-fields  of  these  z-directed  dipoles  are  given  by 


En  =  J 


kn  Idl  e 


-jkr 


4?t 


sinO 


(31) 


.  k  Kdl 

J  4  it 


sinO 


(32) 


With  these  equations  and  the  diffraction  coefficient  developed  in  the 
previous  section  the  computer  code  for  the  radiation  from  the  sources 
shown  in  Figure  4  has  been  developed  in  a  straightforward  manner.  The 
total  field  E*  due  to  the  electric  dipole  at  yn  is  the  sum  of  the  direct 
(or  incident)  field  from  the  dipole  given  by  (31)  with  r  =  rn,  9  =  90° 
and  the  field  diffracted  from  the  cylinder 


E 


d  -  Jkn 

ft  4n 


Idl 


Ds(4>-4n)e 


jkp  cos(4>-<i>  )  "dkrn 
n  n  e 


(33) 


in  which  pn  is  the  distance  between  the  dipole  at  yn  and  the  axis  of  the 
cylinder,  <j>n  is  the  angle  p„  makes  with  the  x-axis,  rn  is  the  distance 
between  the  dipole  at  yn  and  the  field  point,  and  DsU-<(>n)  is  calculated 
for  e  =  ed  =  90°.  Combining  (31)  and  (33) 


rt  .  kn  Idl 

t0  ■  J  4tt 


[1  + 


Ds(<0-fn)e 


jkPncos(<t-<|)n) 


] 


(34) 


r  -  The  total  far-field  due  to  all  the  z-directed  electric  current  dipoles  in 

[  _  the  aperture  is  calculated  by  summing  (34)  on  n  to  obtain  the  contributions 

from  the  2b/w  dipoles  involved. 


IV.  CYLINDER  BLOCKAGE  EFFECTS  ON  THE  LAMPS  ANTENNA 


Let  us  now  proceed  to  use  the  previously  presented  method  to  predict 
the  effects  of  mast  blockage  on  the  LAMPS  antenna.  The  problem  geometry 
is  shown  in  Figure  6.  The  LAMPS  antenna  itself  is  an  offset  fed  para¬ 
bolic  reflector  antenna  with  the  following  estimated  parameters: 

Diameter:  34  inches 

Gain:  28  dB  @  4.6  GHz  (Midband) 

Polarization:  Linear  Vertical 
Beamwidth:  5°  x  5°. 

In  order  to  simplify  our  model  of  the  LAMPS  antenna  we  will  assume 
that  the  effects  of  aperture  blockage  due  the  feed  and  its  supports  can 
be  neglected.  This  assumption  will  have  little  effect  on  our  calculations 
when  cylinder  blockage  is  a  significant  factor,  since  the  pattern  degradation 
caused  by  the  cylinder  will  be  much  greater  than  that  caused  by  aperture 
blockage.  However,  patterns  calculated  with  the  blocking  cylinder  absent 
will  have  far  out  sidelobes  at  an  unrealistically  low  level  [7],  This  should 
be  kept  in  mind  when  evaluating  the  effects  of  the  cylinder  blockage  on 
the  pattern. 

For  lack  of  more  detailed  information  we  will  also  assume  that  the 
aperture  distribution  is  radially  symmetric,  and  of  the  type 


where  b  is  the  aperture  radius  and  here  r  is  measured  radially  in  the  plane 
of  the  aperture  from  its  center  (see  Figure  5).  This  assumption  yields,  for 
the  LAMPS  aperture,  at  4,6  GHz,  the  following  Table  [6]: 


p 

Gain 

1/2  Power  Beamwidth 

0 

32.33 

4.42° 

1 

31 .08 

5.5° 

2 

29.81 

6.37° 

3 

28.76 

7.15° 

Since  the  gain  of  the  antenna  is  more  easily  affected  by  spillover, 
blockage,  losses,  etc.,  than  the  pattern  shape,  we  will  choose  p  =  1 , 
because  this  distribution  produces  the  beamwidth  which  is  closest  to  the 
given  value  of  5°. 

We  can  now  proceed  to  use  (29)  and  (30)  to  evaluate  the  relative 
dipole  magnitudes  Idl  and  Kdl .  Since  the  antenna  is  vertically  polarized, 
Kdl  £  0.  Referring  to  the  geometry  of  Figure  5,  we  have  from  Equation  (29), 
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(35) 


Idl 


2  2  \ 

V  +2  )ejkz  COS 

b2  J 


dy  dz 


2  2  2 

since  r  =y  +z  .  The  factor  C  is  an  unknown  constant  which  need  not  concern 
us  since  we  are  interested  in  relative  pattern  levels  rather  than  absolute 
field  strengths.  Since  the  cylinder  and  the  aperture  dipoles  are  vertical, 
h  =  90°.  If  we  make  the  approximation  that  r^  =  y2  +  z2  (since  w  <<  b), 
and  note  that  z^  =  b2  -  y2,  we  obtain 


Idl  =  C  2w  b2-y2 

4  Jn 


2f2  [ ( b2-y2 ) 3/2+3y2 ( b2-y2 ) 1 /2 ] 


(36) 


for  the  dipole  strength  as  a  function  of  strip  width  w  and  dipole  position 
yn.  Note  that  the  first  term  in  (35)  is  the  dipole  strength  for  a  uniform 
distribution;  the  second  term  is  due  to  the  assumed 


taper.  The  result  obtained  in  (36)  was  used  in  the  previously  described 
computer  program.  Convergence  tests  were  made,  and  it  was  found  that  40 
strips  gave  accurate  results  for  the  13.23a  diameter  LAMPS  antenna.  As 
will  be  seen  the  half  power  beamwidth,  the  1st  side-lobe  level  of  -24.6  dB, 
and  the  1st  null  position  of  7.1°  of  the  unblocked  antenna  pattern  are 
quite  accurately  reproduced  [6]. 

The  calculated  radiation  patterns  in  Figures  7  through  14  will  be 
discussed  next.  In  Figures  7  and  8  the  variation  of  the  blockage  effects 
with  the  mast-to-LAMPS  antenna  separation  is  shown.  The  separation  pp  is 
varied  within  the  given  range  of  2  ft.  to  10  ft.  Note  that  despite  tne 
5  to  1  variation  in  separation  the  on-axis  blockage  varies  less  than 
0.4  dB.  This  occurs  because  we  are  within  the  near  field  of  the  antenna 
throughout  this  range  of  p0  variation.  Indeed,  using  the  D2/>  far  field 
criteria,  the  near  field  for  the  LAMPS  antenna  extends  to  37.5  feet. 

Thus  for  the  range  of  p0  values  indicated  in  the  Figures  7  and  8  the 
antenna  beam  has  not  started  to  spread  and  the  1/r  spherical  wave  field 
behavior  is  not  yet  evident.  Thus  any  attempt  to  reduce  the  on  axis 
blockage  by  increasing  the  antenna-cylinder  spacing  p0  will  have  little 
effect  for  spacings  of  less  than  30-40  feet.  The  blockage  effect  for  on- 
axis  separations  of  D2/a  =  37.5  feet  and  2D2/a  =  75  feet  are  shown  in 
Figure  9.  The  scattered  field  from  the  cylinder  is  known  to  decrease  by 
3  dB  when  the  separation  between  the  antenna  and  cylinder  is  doubled  in 
the  far-zone  of  the  antenna.  This  decrease  is  apparent  in  the  region 
*>40",  where  the  field  directly  radiated  from  the  antenna  is  insignificant. 
The  small  irregularities  in  the  two  patterns  are  caused  by  the  large 
separation  of  the  cylinder  and  aperture  with  a  resulting  rapid  phase 
variation  with  aspect  $  of  the  second  term  within  the  brackets  of  (34). 
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In  Figures  10  and  11  the  variation  of  blockage  effects  with  cylinder 
radius  is  shown.  The  larger  cylinders  have  greater  on-axis  blockage  due 
to  their  greater  shadowing  effect.  The  exact  on-axis  blockage  loss  for 
the  curves  of  figures  10  and  11  are  as  indicated  in  the  figure  titles. 

In  Figures  12,  13  and  14  the  effects  of  varying  the  angular  position 
of  the  blocking  cylinder  are  illustrated.  In  Figure  12  patterns  are  cal¬ 
culated  for  <P0  values  of  0°,  10°,  and  30°.  For  <J>0  =  10°  the  cylinder  is 
near  the  edge  of  the  near  field  "beam11  and  the  blockage  is  greatly  de¬ 
creased  from  the  4>0  =  0°.  For  <J>0  =  30°  the  cylinder  is  well  out  of  the 
active  near  field  region,  and  the  scattering  effects  of  the  cylinder  are 
negligible  except  for  the  far-out  sidelobe  levels. 

Figures  13  and  14  confirm  that  a  knowledge  of  the  far  field  pattern 
of  the  LAMPS  antenna  does  not  provide  insight  into  near  field  blockage 
effects.  For  example,  the  far  field  pattern  has  a  null  at  <J>  =  7°  and  a 
first  sidelobe  at  4>  =  9°.  Thus,  were  the  cylinder  in  the  far  field  of 
the  antenna,  locating  the  cylinder  at  f>0  =  7°  should  produce  less  effect 
than  with  the  cylinder  at  4>o  =  9°.  The  results  shown  in  Figures  13  and 
14  indicate  that  the  opposite  results,  and  that  the  far  field  distrioution 
has  not  yet  been  formed  for  a  separation  of  only  6  feet. 

At  this  point,  as  a  reminder,  we  repeat  that  the  unperturbed  patterns 
shown  in  Figures  7  through  14  were  calculated  without  including  blockage 
effects  of  the  feed  system;  thus  the  sidelobe  levels  are  unrealistically 
low.  For  a  more  valid  comparison  one  should  compare  the  calculated  side¬ 
lobe  levels  of  the  perturbed  (by  the  cylinder)  patterns  with  measured 
patterns  of  the  LAMPS  antenna.  However,  the  figures  contained  herein 
do  provide  a  good  estimate  of  the  effects  of  the  cylinder,  since  for  the 
range  of  p0,  <t>o>  and  a  values  considered  the  blockage  effects  of  the 
cylinder  are  usually  much  more  significant  than  those  of  the  feed  system. 
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Figure  6 
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SIDE  VIEW 


LAMPS  antenna  with  a  blocking  cylinder. 
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pattern  of  the  LAMPS  antenna  with  a  blocking  cylinder 


blocki 


V. 


CONCLUSIONS 


An  expression  for  the  cylindrical  diffraction  coefficient  suitable  for 
use  with  circular  cylinders  has  been  developed.  The  diffraction  coefficient 
is  in  the  GTD  format,  and  can  be  incorporated  into  a  general,  ray-optical 
analysis  of  a  radiating  system.  Using  this  diffraction  coefficient  a 
computer  code  has  been  written  which  is  capable  of  calculating  antenna 
patterns  for  an  aperture  antenna  in  the  vicinity  of  a  circular  conducting 
cylinder.  At  present  the  cylinder  must  be  parallel  to  the  plane  of  the 
aperture,  and  the  pattern  cut  must  be  in  the  plane  perpendicular  to  the 
cylinder  axis.  However,  the  theory  needed  for  an  extension  of  the  code 
to  the  more  general  case  is  included  herein.  Recently  these  diffraction 
coefficients  have  been  used  to  calculate  the  aperture  blockage  in  the  case 
of  a  reflector  antenna  with  cylindrical  feed  supports  [7]. 

Using  this  computer  code  the  blockage  effects  of  a  3"  to  6"  radius 
mast  located  2  to  10  feet  from  the  34"  diameter  LAMPS  antenna  were  cal¬ 
culated  at  4.6  GHz.  It  was  found  that  on-axis  blockage,  with  the  cylinder 
directly  in  front  of  the  antenna  ranged  from  -3.67  dB  for  a  3"  radius 
cylinder  to  -7.96  dB  for  a  6"  radius  cylinder  with  the  distance  held 
constant  at  6'.  The  on-axis  blockage  was  affected  very  little  by  changing 
the  antenna  mast  separation  PQ,  staying  with  ±0.2  dB  of  -4.8  dB  as  PQ 
varied  from  2'  to  10',  with  a  4"  cylinder  again  located  directly  in  front 
of  the  antenna.  This  insensitivity  to  mast-antenna  separation  is  due  to 
the  fact  that  the  mast  is  located  well  within  the  near  field  region  of  the 
LAMPS  antenna  for  separations  less  than  30-40  feet.  The  on-axis  blockage 
was  found  to  decrease  in  the  appropriate  manner  as  p0  was  increased  from 
37.5  to  75  feet. 

With  p0  held  at  6  feet  the  on-axis  blockage  was  reduced  significantly 
as  the  mast  was  rotated  from  in  front  of  the  antenna.  As  i>0  changed  from 
0°  to  10°  to  30°  the  on-axis  blockage  was  reduced  from  -5.01  dB  to  -1.416  dB 
to  0.0  dB.  In  general,  the  amount  of  angular  displacement  of  the  mast  away 
from  on  axis  (<}>0  =  0°)  required  to  reduce  the  blockage  below  a  given  level 
is  reduced  by  increasing  P0  the  antenna-mast  separation. 
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APPENDIX 

THE  FIELDS  OF  AXIAL  DIPOLES  IN  THE  PRESENCE  OF  A  CIRCULAR  CYLINDER 


The  geometry  and  coordinates  are  shown  in  Figure  1.  The  Green's 
function  method  will  be  used  to  find  the  radiated  fields.  Due  to  the 
nature  of  the  problem,  we  need  only  find  Ez,  the  z  component  of  the 
electric  field,  for  the  electric  dipole  case,  and  Hz  for  the  magnetic 
dipole  case.  As  in  the  text,  the  cylinder  is  perfectly  conducting.  One 
these  components  are  found  then 


E0  =  -Ez/sin  e 

(Al) 

E  =  nH  /sin  6 
<P  z 

(A2) 

The  scalar  Green’s  functions  of  interest 
geneous  partial  differential  equation 

here  satisfy  the  inhomo- 

(v2  +  k2)  G(r,r' )  =  -6(7-7') 

(A3) 

together  with  the  radiation  condition  at  infinity  and  certain  boundary 
conditions  on  S.  The  Green's  function  G-|(r,r')  =  0  when  r  is  on  the 
cylindrical  surface,  whereas 

3G_  (r,r' ) 

3TT  =  0 

when  r  is  on  the  cylindrical  surface.  Here  n'  is  taken  in  the  direction 
normal  to  the  cylindrical  surface.  For  the  circular  cylinder 

3  3_ 

in"1"  3p  ’ 

r  is  the  position  vector  for  the  field  point  (p,<t>,z),  and  r'  is  the 
position  vector  for  the  source  point  (p  '  ,<f> 1  ,z ' ) . 


It  can  be  shown  that  the  solution  to  the  boundary  value  problem  for 
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in  which 

i  e~ jkR 

G  (R)  =  rr-  p -  is  the  free-space  Greer  's  function, 

R  =  ; r - r 1 S  is  the  distance  between  the  field  and  source  points, 

►  =J  kZ-t?  , 

and  the  subscripts  ■  ,  *  mean  choose  the  value  of  r  or  p*  whichever 
is  smaller  or  larger,  respectively . 

In  the  far-zone  one  may  replace  Hj  (•  ■  )  by  its  large  argument 
asymptotic  approximation;  furthermore  (z-z1)  is  large,  so  the  above 
integral  in  Ge  may  be  evaluated  asymptotically  to  obtain 


G  ( r, r 1  ) 
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In  this  discussion  the  source  is  the  electric  current  density 


Jv(r* )  =  pe  «(r'-p0)z  , 


(A6) 


where  the  position  vector  p0  defines  the  source  location  (p0,o,o)  shown 
in  Figure  1 . 


that 


For  such  a  z-directed  source  it  follows  from  Maxwell's  equations 


[v2  +  k2]E  =  jcou 


1  + 
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k2  3z2 


]v 


(A7 ) 


Using  Green's  second  identity  or  comparing  the  above  equation  with  (A3), 
it  follows  that 


Ez(r)  =~jwu 


1  + 


1  o 


k2  3Z2 


Gl(|"'‘,o>Pe 
In  the  far-zone  this  simplifies  to 


(A8) 


Ez(r)  =-jkn  sin20  G1(r,PQ)pe, 


(A9) 


where  wu  has  been  replaced  by  kn.  Thus  from  (Al),  (A4),  (A5)  and  (A9) 
and  noting  that 


R  a  r  -  po  sinO  cos* 


with  r  >>  pq  in  the  far-zone, 
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which  is  the  same  as  (4)  in  the  text. 


The  magnetic  current  moment  (dipole)  pm  gives  rise  to  a  .‘-directed 
electric  field.  The  derivation  of  the  expression  for  D  closely  oarallel 
the  procedingdevelopment.  One  starts  with  G2(r,r')  which  has  the  same 
form  as  G](r,r')  except  that  the  ratio 
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where  the  prime  denotes  di fferentiation  with  respect  to  the  argument, 
Then  following  the  steps  taken  in  the  preceding  case  and  employing 
duality,  we  find  that 
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and  so  from  (A2) 
r  -  'Ik  n 
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in  the  far-zone.  The  first  term  in  the  above  equation  is  the  field 
directly  radiated  from  the  magnetic  dipole  (the  incident  field),  and 
the  second  term  is  the  field  diffracted  by  the  cylinder.  The  second 
term  is  used  to  find  the  hard  diffraction  coefficient  and  so  to  define 
in  (22). 


